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(B) AMENDMENTS TO THE CLAIMS 

1 . (Currently amended) A method for seismic data processing, comprising: 

constructing explicit depth extrapolation operators with variable operator lengths 
depending on maximum dip angle, accuracy condition, and wavenumber; aad 
constructing operator tables using the explicit depth extrapolation operators having the 
smallest operator length satisfying the accuracy condition at the maximum dip angle for 
each of a plurality of wavenumbers; anfl 

performing depth migration using the _explicit depth extrapolation operator in a sub- 
domain at a depth having the smallest operator length from the operator tables for the 
highest wavenumber in the sub-domain at the depth. 

2. (Canceled) 

3. (Original): The method of claim 1, wherein the step of constructing operator tables 
comprises: 

selecting a maximum dip angle; 
selecting a type of operator, 

selecting an accuracy condition for the type of operator; 
selecting a plurality of wavenumbers; and 

performing the following steps for each of the plurality of wavenumbers: 
selecting a plurality of operator lengths for the selected wavenumber; and 
performing the following steps for each of the plurality of operator lengths: 

detemiining if the selected type of operator with the selected operator length 
satisfies the selected accuracy condition at the selected maximum dip angle at the 
selected wavenumber; 

determining if the operator length is the smallest operator length satisfying the 
accuracy condition at the selected maximum dip angle for the selected 
wavenumber, and 
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storing the operator length in an operator table if the operator length is the 
smallest operator length satisfying the accuracy condition at the maximum dip 
angle for the selected wavenumber. 

4. (Previously amended): The method of claim 3, wherein the steps performed for each 
wavenumber comprise: 

selecting a plurality of pairs of first operator half-lengths and second operator half- 
lengths for the selected wavenumber, and 

performing the following steps for each of the plurality of pairs of first operator half- 
lengths and second operator half-lengths: 

determining if the selected type of operator with the selected pair of first operator 
half-length and second operator half-length satisfies the selected accuracy condition 
at the selected maximum dip angle at the selected wavenumber, 
determining if the first operator half-length is the smallest first operator half-length 
satisfying the accuracy condition for the selected wavenumber and the second 
operator half-length is the smallest second operator half-length satisfying the 
accuracy condition for the selected wavenumber; and 

storing the pair of first operator half-length and second operator half-length in an 
operator table if the first operator half-length is the smallest first operator half-length 
satisfying the accuracy condition for the selected wavenumber and the second 
operator half-length is the smallest second operator half-length satisfying the 
accuracy condition at the selected maximum dip angle for the selected wavenumber. 

5. (Previously amended): The method of claim 1, wherein the step of constructing operator 
tables comprises: 

selecting a maximum dip angle; 
selecting a type of operator; 

selecting an accuracy condition for the type of operator; 

determining preliminary operator tables with variable operator lengths; 

selecting a plurality of wavenumbers; and 
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performing the following steps for each of the plurality of wavenumbers: 

selecting a smallest operator length in the preliminary operator tables satisfying the 
accuracy condition for the selected wavenumber; and 
storing the smallest operator length in an operator table. 

6. (Previously amended): The method of claim 5, wherein the steps performed for each 
wavenumber comprise: 

selecting a pair of smallest first operator half-length and smallest second operator half- 
length in the preliminary operator tables satisfying the accuracy condition for the selected 
wavenumber; and 

storing the pair of smallest first operator half-length and smallest second operator half- 
length in an operator table. 

7, (Previously amended): The method of claim 2 9 wherein the step of performing depth 
migration comprises: 

selecting a seismic data set; 

selecting a plurality of depths in the seismic data set; 
selecting sub-domains for each of the plurality of depths; 
determining a velocity model for the seismic data set and the plurality of depths; 
determining a lowest velocity for each sub-domain from the velocity model; 
selecting an operator table with variable operator lengths; 
selecting a plurality of frequencies in the selected sub-domain; and 
performing the following steps for each of the plurality of frequencies: 
performing the following steps for each of the plurality of depths: 

performing the following steps for each of the subnlomains at the depth: 

calculating a highest wavenumber in the sub-domain from the frequency and 

the lowest velocity for the sub-domain; 

selecting a maximum operator length for the calculated highest wavenumber 
from the operator table; 
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applying operators from the operator table with operator lengths no greater 
than the maximum operator length to the selected sub-domain. 

8. (Previously amended): The method of claim 7. wherein the operator table with variable 
operator lengths comprises explicit extrapolation operators with variable operator lengths. 

9. (Previously amended): The method of claim 7, wherein the step of selecting an operator 
table with variable operator lengths further comprises: 

interpolating the operator table with variable operator lengths. 

10. (Original): The method of claim 7, wherein each operator length comprises a pair of first 
operator half-length and second operator half-length* 

11. (Original); The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 

P(x i ,y j9 a>,z + Az) = 

where P(x b y Jt &, z+Az) is a seismic wavefield at lateral location * f » iAx, yj = jAy f and depth 
z+Az; o) is angular frequency; Pfa yj, co, z) is a seismic wavefield in a sub-domain D at depth z; 
Ax, Ay, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
L(k Qf max (D l z)) and M(k^ mQX (D t z)) are operator half-lengths in the x- and y- coordinate directions, 
respectively, for the highest wavenumber kJ^fD, z) in the sub-domain D at depth z\ Wfa y Jt 
kjxu y Jt z) t Az) is an explicit depth extrapolation operator; and k^fa, y Jt z) is a local wavenumber. 
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12. (Original): The method of claim 11, wherein the wavenumber kjx^ y Jf z) is given by the 
formula: 

J c(x n yj,z) 
where c(x u yj> z) is a local propagation velocity, 

13. (Original): The method of claim 11, wherein the highest wavenumber kv max (D> z) for a 
frequency co for a sub-domain D at depth z is given by the formula; 

where c min (D, z) is the lowest velocity for the sub-domain D at depth z, 

14. (Original): The method of claim 3, wherein the step of step of selecting a type of operator 
comprises applying the following formula: 

P(x i9 a> y z + Az) = J] w (*i> k *( x t> z )> p &i 

where P(xu <o, z+Az) is a seismic wavefield at lateral location x$ m i-Ax and depth z+Az, a> is 
angular frequency, Pfa <o t z) is a seismic wavefield in a sub-domain D at depth z 9 Ax and Az are 
step lengths in the and 2-coordinate directions, respectively, LfkJ^fD. z)) is an operator half- 
length in the x-coordinate direction for the highest wavenumber kco max (D, z) in the sub-domain D 
at depth z; W(x if k^(x it z), Az) is an explicit depth extrapolation operator, and k^fa, z) is a local 
wavenumber. 
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15. (Original): The method of claim 14, wherein the wavenumber kjfa z) is given by the 
formula: 



where c(x b z) is a local propagation velocity. 

16. (Original): The method of claim 3, wherein the step of selecting a type of operator comprises 
applying the following formula: 



where P(x if y^ a>, z+Az) is a seismic wavefield at lateral location x t « iAx t yj = y'-Ay, and depth 
o) is angular frequency; Pfa y^ a>, 3) is a seismic wavefield in a sub-domain D at depth z\ 
Ax , Ay, and Az are step lengths in the x-, y- and z-coordinate directions, respectively; 
L(kc max (D r z)) is an operator length for the highest wavenumber ka> max (D> z) in the sub-domain D 
at depth z; Wi(k^(xu yj, z), Az) is an explicit depth extrapolation operator expressed in cylindrical 
coordinates; k^fe yj> z) is a local wavenumber; and htfe y jt w f z) ^ are auxiliary fields resulting 
from applying recursive 2D Chebyshev filters to the seismic wavefield P(xu yj, co> z). 



k„(x f ,z) = 



c(x n z)' 



P(x i ,y J ,o>z + Az) = }L W f {k m (x n y Jt z) 9 6z)k 9 {x n y t 



a>,z) t 
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